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EXECUTIVE  SUMMARY 

The  Fram  Strait  Acoustic  Monitoring  for  Ocean  Climate  Studies  Workshop  was  held  at 
the  NOAA  Environmental  Technology  Laboratory  in  Boulder,  Colorado,  on  September  25-26, 
1996.  The  workshop  was  convened  to  discuss  and  develop  a  prioritized  scientific  plan  focused 
on  a  systems  approach  to  the  oceanographic  investigation  of  the  Fram  Strait.  Major  topics 
discussed  during  the  workshop  included  the  following: 

•  monitoring  of  heat  and  water  mass  influx  through  the  Fram  Strait 
as  part  of  a  general  problem  of  Arctic  warming 

•  acoustic  tomography  monitoring  of  the  Fram  Strait 

•  monitoring  of  flux  into  the  Arctic 

•  planning  of  a  pilot  experiment  and  feasibility  test 

The  following  issues  and  questions  emerged  from  the  workshop  deliberations: 

(1)  The  Arctic  Ocean  is  a  system  of  high  climate  sensitivity.  An  enhanced  greenhouse 
warming  in  the  Arctic  relative  to  lower  latitudes  by  a  factor  of  2-4  will  have  a  drastic  impact  on 
ice-pack  shrinking  and  the  global  climate. 

(2)  The  heat  flux  exchange  between  the  Arctic  and  the  North  Atlantic  at  the  Fram  Strait 
plays  a  significant  role  in  the  global  climate  studies  for  ATOC-Arctic.  Acoustic  monitoring  of 
the  heat  flux  at  the  Fram  Strait  is  proposed,  and  issues  related  to  the  acoustic  propagation  in  this 
extremely  complex  environment  are  investigated. 

(3)  Different  methods  of  acoustic  temperature  and  current  measurements  will  be 
considered,  including 

•  arrival-time  variation  for  temperature  and  along-path  current  measurement 

•  parametric  array  for  selective  mode  excitation  to  increase  the  sensitivity  of  the  method 

•  scintillation  technique  for  across-path  current  measurement 

•  horizontal  refraction  tomography  for  current  and  mesoscale  structure  identification 

Acknowledgments.  We  would  like  to  thank  all  the  participants  of  the  workshop  for  their  active 
involvement  during  the  meeting.  In  addition,  we  appreciate  the  assistance  of  Sharon  Kirby-Cole 
and  Karen  Martin  during  the  workshop  and  in  report  preparation.  Funding  for  the  workshop  was 
cosponsored  by  the  NOAA  Environmental  Technology  Laboratory  (USA)  and  the  Nansen 
Environmental  and  Remote  Sensing  Center  (Norway).  This  project  was  partially  supported 
NATO  Linkage  Grant  (E.NVIR.L  G  960352)  528  (96)  LVdC). 


1.  INTRODUCTION 

In  a  meeting  in  Rome  in  December  1995,  the  Intergovernmental  Panel  on  Climate  Change 
(IPCC)  concluded  that  the  observed  global  warming  of  0.4°-0.6°C  in  the  last  century  is  partly 
caused  by  human  activity  (Hasselmann,  1993;  IPCC,  1996a-c).  The  increasing  emission  of 
greenhouse  gases  into  the  atmosphere  is  predicted  to  cause  an  average  global  warming  in  the  next 
century  of  approximately  3°C  (Cubasch  et  al.,  1995;  Mitchell  et  al.,  1995)  with  enhanced  effects 
in  the  Arctic  region.  Due  to  the  lack  of  an  integrated  monitoring  system  of  oceanographic 
parameters  sensitive  to  climate  change,  apart  from  sea-ice  monitoring  from  satellites 
(Johannessen  et  al.,  1995a,  1996),  there  are  no  adequate  observations  for  estimating  the  large- 
scale  ocean  temperature  variability,  which  can  be  considered  as  a  noise  background  where  the 
greenhouse  gas-induced  warming  signal  must  be  detected.  Today,  we  have  to  rely  on  model 
estimates,  based  on  poor  observational  records,  for  both  the  global  warming  signal  and  the 
ambient  background  climate  noise. 

Climate  modeling  studies  have  suggested  that  the  Arctic  region  will  show  the  strongest 
atmospheric  warming  in  response  to  increased  greenhouse  effect  (Manabe  et  al.,  1991;  Cattle  and 
Crossley,  1995).  This  may  cause  a  melting  of  the  sea  ice  (Semtner,  1987;  Johannessen  et  al., 
1995a,b)  and  a  warming  of  the  Atlantic  water  being  advected  through  the  Fram  Strait  and  the 
St.  Anna  Trough,  causing  an  increase  of  the  internal  temperature  of  the  Arctic  Ocean 
(Macdonald,  1996).  Therefore,  the  Fram  Strait,  with  inflow  and  outflow  from  the  Arctic  Ocean, 
is  a  key  area  to  study  the  integrated  effect  of  global  warming  both  on  the  Nordic  Seas  and  on  the 
Arctic  Ocean  itself  (Fig.  1). 

The  Fram  Strait  is  the  region  where  almost  all  heat  and  water  exchange  between  the 
Arctic  Ocean  and  the  Atlantic  Ocean  occurs,  apart  from  the  Atlantic  inflow  through  the  St.  Anna 
Trough.  To  a  first-order  approximation,  the  Strait  acts  both  as  an  entry  and  an  exit  port  for 
the  Arctic  Ocean.  The  general  large-scale  ocean  circulation  in  this  region  is  dominated  by  the 
shallow,  southward-flowing,  cold,  low-salinity  East  Greenland  current  (EGC)  that  exports  ice 
and  polar  water  out  of  the  Arctic  Ocean,  and  the  northward-flowing,  warm,  saline  Atlantic  water 
in  the  West  Spitsbergen  current  (WSC),  shown  in  Fig.  2.  Heat  inflow  through  the  Fram  Strait  in 
proportion  to  the  total  is  about  80%  (see  Table  1). 

The  water  mass  exchanges  through  the  Fram  Strait  are  essential  boundary  conditions 
for  large-scale  ocean  circulation  and  global  climate  models.  In  recognition  of  this  fact,  the 
observational  focus  in  the  Fram  Strait  region  has  been  on  measuring  the  fluxes  of  mass  and  heat 
through  the  Strait  in  order  to  better  understand  the  difficult-to-observe  water  mass  structure,  due 
to  the  ice,  of  the  Arctic  Ocean  (Aagaard  and  Carmack,  1994).  Often,  estimates  have  been  made 
combining  current  data  and  hydrographic  data  obtained  at  different  times  and  locations.  This  has 
resulted  in  a  large  spread  of  the  estimated  values.  Previous  estimates  of  volume  transport  of  the 
WSC  across  79°  N  vary  from  1.9x10^  to  8x10^  m^  s"'  (Simonsen  and  Haugan,  1996).  This  large 
variation  in  water  mass  exchange  introduces  errors  in  ocean  circulation  models. 
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Figure  1.  The  Fram  Strait  is  a  key  area  to  study  the  integrated  effect  of  global  warming. 
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Figure  2.  Ocean  current  structure. 


Table  1.  Heat  Flux  Exchange  Rate  through  the  Fram  Strait^ 


Volume  Transport 

Heat 

Transport 

Mean  Temperature 

(SJ 

(10' 

'kcalS') 

CO 

East  Greenland  current 

Polar  water 

-1.8 

2.0 

-1.2 

Atlantic  water 

-5.3 

-3.2 

0.5 

Ice 

-0.1 

8.0 

West  Spitsbergen  current 

7.1 

16.3 

2.2 

Total  inflow 

9.4 

1.8 

Total  outflow 

-9.4 

— 

0.1 

Total  advective  heat 

29.7 

Total  advective  heat  loss 

-3.9 

Net  exchange  0.0  25.8  — 

Heat  inflow  of  the  Fram  Strait  in  proportion  to  the  total:  80% 
^Aagaard  and  Greisman  (1975). 

The  main  problem  in  the  investigation  of  straits  in  general  is  the  lack  of  synoptic  and 
high-temporal  resolution  hydrographic  and  current  measurements  to  estimate  the  volume  and 
heat  flux  through  straits.  The  objective  of  the  Fram  Strait  monitoring  project  is  to  measure 
synoptic  temperature  and  inflow  and  outflow  in  order  to  improve  estimates  of  volume  and  heat 
flux  through  the  Fram  Strait.  Acoustic  techniques  can  play  an  important  role  in  this  monitoring 
because  of  the  spatial  averaging  that  is  inherent  in  the  methods. 

Both  passive  and  active  methods  in  underwater  acoustics  have  become  important  in 
detection,  exploration,  and  communication  systems  for  use  in  oceanographic  investigations  and 
in  the  oil  industry.  A  promising  method,  and  probably  the  only  one  to  measure  gyre-scale 
average  ocean  temperature,  is  to  use  acoustic  monitoring  based  on  the  concept  of  Acoustic 
Thermometry  of  Ocean  Climate  (ATOC;  Munk  and  Forbes,  1989).  This  concept  is  an  extension 
of  acoustic  tomography,  also  pioneered  by  Walter  Munk  (Munk,  1994;  Munk  et  al.,  1995). 
Acoustic  tomography  has  been  successfully  used  to  study,  for  example,  ocean  mesoscale 
(Comuelle  et  al.,  1985),  a  small-scale  ocean  structure  (Comuelle  and  Howe,  1987),  and  the  deep 
water  convection  in  the  Greenland  Sea  (Worcester  et  al.,  1993),  and  in  other  experiments 
(Comuelle  et  al.,  1985;  Comuelle  and  Howe,  1987).  The  acoustic  thermometry  concept  was 
tested  successfully  in  the  Heard  Island  Feasibility  Test,  conducted  in  early  1991  (Munk  et  al., 
1994;  Heard  Island  Papers,  1994).  During  a  pilot  Transarctic  Acoustic  Propagation  (TAP) 
experiment  conducted  by  MikhaJevsky  et  al.  (1994,  1995)  in  the  spring  of  1994,  long-range 
propagation  (2700  km)  through  the  Arctic  Basin  was  performed.  The  TAP  experiment 
demonstrated  that  transmission  across  the  Arctic  Basin  is  possible  at  a  source  frequency  of  20  Hz. 
Some  data  about  the  ice  cover  effect  on  the  sound  propagation  were  obtained  by  Jin  et  al.  (1993, 
1994). 


A  series  of  experiments  of  steadily  increasing  ranges  from  25  to  1275  km  have  explored 
the  feasibility  of  using  reciprocal  acoustic  transmission  to  measure  large-scale  ocean  current. 
In  the  1987  Reciprocal  Tomography  Experiment,  located  in  the  central  North  Pacific  Ocean,  the 
eight  tidal  current  components,  amplitudes,  and  phases  calculated  by  acoustic  methods  were  in 
excellent  agreement  with  barotropic  tidal  current  computed  from  a  current  meter  mooring 
(Luther  et  al.,  1991;  Munk  et  al.,  1995);  further  excellent  agreement  has  been  found  with  the 
Ocean  Topography  Experiment  (TOPEX)  altimetry  data  (Dushaw  et  al.,  1997).  A  first  attempt  to 
study  the  possibility  of  using  two-dimensional  acoustic  propagation  to  monitor  a  stratified  shear 
flow  was  made  by  Farmer  (Farmer  and  Di  lorio,  1996).  The  experiment  was  designed  to  study 
the  fine  structure  of  the  turbulent  current  in  a  3-km-long,  1-km-wide,  30-m-deep  channel  using 
high-frequency  sound  (67  kHz).  Based  on  this  investigation,  they  concluded  that  two- 
dimensional  acoustic  arrays  can  provide  an  effective  tool  for  studying  the  stratified  shear  flows  in 
coastal  regions,  especially  for  determining  the  component  of  the  flow  perpendicular  to  the 
acoustic  path.  The  philosophy  of  using  acoustic  arrays  (sources  and  receivers)  to  measure  the 
current  structure  can  be  transferred  to  similar  concepts  for  a  much  larger  strait,  such  as  the  Fram 
Strait,  using  frequencies  below  200-500  Hz  and  larger  array  configurations  (see  Fig.  1).  Effects 
of  ice  cover  on  the  acoustic  signal  propagation  (Jin  et  al.,  1993,  1994)  have  to  be  taken  into 
account. 

Based  on  the  above  scientific  rationale,  the  Scientific  Committee  on  Oceanic  Research 
Working  Group  (SCOR  WG)  96  made  the  following  recommendations: 

"The  SCOR  WG  96  recognizes  the  particular  sensitivity  of  high-latitude  regions  to 
global  wanning  and  identifies  the  desirability  of  an  early  implementation  of  an  Arctic 
acoustic  monitoring  program." 

A  novel  concept  for  an  Arctic  Ocean  Climate  Monitoring  System  is  outlined  in  Fig.  3 
(Mikhalevsky  et  al.,  1994,  1995).  It  represents  an  innovative,  unique  system  combining  acoustic 
remote  sensing  of  the  ocean  and  remote  sensing  of  ice  from  space  (Johannessen  et  al.,  1995a), 
including  modeling  and  data  assimilation  (Evensen,  1994a,b).  However,  in  this  project  the  focus 
is  on  the  feasibility  of  using  acoustic  techniques  for  monitoring  volume  and  heat  flux  through  the 
Fram  Strait,  as  indicated  by  the  shading  of  the  Strait  component  in  Fig.  3. 

Summing  up,  we  conclude  the  following: 

(1 )  The  Arctic  Ocean  is  a  system  of  high  climate  sensitivity.  An  enhanced  greenhouse 
warming  in  the  Arctic  relative  to  lower  latitudes  by  a  factor  of  2-i  will  have  a  drastic  impact  on 
ice-pack  shrinking  and  the  global  climate. 

(2)  The  heat  flux  exchange  between  the  Arctic  and  the  North  Atlantic  at  the  Fram  Strait 
plays  a  significant  role  in  the  global  climate  studies  for  ATOC-Arctic.  Acoustic  monitoring  of 
the  heat  flux  at  the  Fram  Strait  is  proposed,  and  issues  related  to  the  acoustic  propagation  in  this 
extremely  complex  environment  are  investigated. 


(3)  Different  methods  of  acoustic  temperature  and  current  measurements  will  be 
considered,  including 

•  arrival-time  variation  for  temperature  and  along-path  current  measurement 

•  parametric  array  for  selective  mode  excitation  to  increase  the  sensitivity  of  the  method 

•  scintillation  technique  for  across-path  current  measurement 

•  horizontal  refraction  tomography  for  current  and  mesoscale  structure  identification 
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Figure  3.  Combined  remote-sensing  system  for  Arctic  climate  monitoring. 
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2.  OCEANOGRAPHIC  MODEL 

2.1  The  Sound  Speed  Field  Structure 

To  model  the  sound  propagation  in  the  Fram  Strait,  a  region  along  79°  N  latitude  could 
be  chosen  where  an  oceanographic  cross  section  was  performed  by  the  RA^  Polarstem 
expedition  in  March  1993.  A  map  of  the  region  with  bathymetry  and  sound  speed  channel  depth 
distributions,  as  well  as  the  suggested  positions  of  the  acoustic  transmitter  (T)  and  receiver  (R), 
is  shown  in  Fig.  4.  Figure  5  shows  the  isolines  of  sound  speed  in  the  vertical  cross  section  of  the 
Strait  along  the  sound  pathway,  based  on  the  experimental  data  obtained  by  RA^  Polarstem,  and 
reflects  the  acoustic  features  of  this  region. 

As  is  seen  in  Fig.  4,  the  larger  spatial  variations  in  sound  speed  on  this  pathway  were 
observed  in  the  upper  layers  of  the  ocean.  The  ocean  inhomogeneities  that  lead  to  changes  in 
sound  speed  can  be  assumed  to  vary  randomly  in  time  and  space.  The  horizontal  spatial  scale 
of  the  sound  speed  inhomogeneity  is  on  the  order  of  30  km,  and  the  sound  speed  can  vary  by 
10  ms"'  at  such  a  spatial  scale.  The  environmental  conditions  at  the  eastern  and  the  western  part 
of  the  pathway  are  different.  The  West  Spitsbergen  current  dominates  at  the  region  adjacent  to 
Spitsbergen,  and  the  East  Greenland  current  dominates  at  the  eastern  part  of  the  acoustic  pathway 
(see  Fig.  2). 

Vertical  profiles  of  sound  speed,  calculated  at  different  points  based  on  the  measurements 
made  by  WW  Polarstem  in  March  1993,  are  presented  in  Fig.  6.  These  profiles  were  considered 
as  initial  steps  to  further  calculations  of  the  sound  propagation  in  a  variable  environment.  Then, 
two  ensembles  of  sound  speed  distributions  were  constructed.  One  was  obtained  by  shifting  the 
sound  speed  field  horizontally  to  simulate  oceanographic  variability  and  to  numerically  estimate 
the  signal  arrival-time  fluctuations.  Another  set  of  profiles  was  constructed  to  model  the  ocean's 
surface-layer  heating. 

2.2  Ocean  Temperature-Increasing  Model 

It  was  assumed  that  the  temperature  changes  take  place  in  the  upper  200-m  layer  of  the 
ocean.  The  corresponding  sound  speed  variations  were  approximated  as  follows: 

—  =  a Ar  , 

c 

,  (1) 

3.09  •  10"^ 

where  Ar  is  an  ocean  water  temperature  variation  in  kelvins,  so  Ar  =  1  K  corresponds  to 
4.51  m  s"'  in  sound  speed  variation  at  c  =  1.46  m  s"'.  The  variations  in  sound  speed  were  added 
to  the  initial  sound  speed  profiles  (SSPs)  in  the  upper  200-m  subsurface  layer: 


THE  SOUND-SPEED  STRUCTURE 


Figure  4.  Regional  bathymetry  and  sound  speed  channel  depth  (Hurdle,  1986),  along  with  a  suggested 
transmitter  (T)  and  receiver  (R)  configuration  for  acoustic  monitoring  of  volume  and  heat  flux. 
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Figure  5.  Sound  speed  isolines  in  the  Fram  Strait  cross  section. 


SOUND  SPEED  PROFILES  AT  DIFFERENT  CROSS  SECTIONS  OF  THE  TRACE 
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Figure  6.  Sound  speed  profiles  across  the  Fram  Strait. 
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6c  =  DC  exp{-z/200}  , 

(2) 
DC  =  4.5ms"V°C  . 

The  modified  SSPs  were  obtained  by  exponential  matching  of  the  profiles  obtained  with  the 
initial  SSPs  in  the  deeper  layer  of  the  ocean. 


3.  ACOUSTIC  PROPAGATION  SIMULATION 

3.1  Ray  Arrival-Time  Pattern 

Inversion  of  acoustical  data  is  a  reasonable  way  of  remote  sensing  large-scale  processes  in 
the  ocean  interior.  Variations  of  travel  time  are  directly  related  to  the  ocean  temperature 
averaged  over  a  corresponding  acoustic  path.  Using  forward-backward  measurements,  one  can 
also  infer  information  about  the  average  velocity  component  along  the  acoustic  path.  To  infer 
the  current  velocity  in  a  direction  perpendicular  to  the  acoustic  path,  the  possible  application  of 
the  scintillation  method  or  horizontal  refraction  method  should  be  investigated.  Consider  first 
the  arrival-time  variation  due  to  ocean  temperature  change  in  the  Fram  Strait  environment 
(Naugolnykh  et  al.,  1997). 

Statistical  ensembles  of  the  sound  speed  fields  were  constructed  from  specific  parts  of 
the  general  hydrographic  data  distribution  presented  in  Fig.  5.  The  space  between  source  and 
receiver  forms  a  stationary  path  of  200-km  length.  Ensembles  of  the  signal  travel  times  were 
obtained  by  shifting  the  entire  sound  speed  pattern  by  10  km.  In  the  course  of  computational 
modeling,  1 1  patterns  of  environmental  conditions  were  considered.  Thus,  the  total 
environmental  shift  in  our  modeling  was  1 10  km,  which  is  essentially  larger  than  the  spatial  scale 
of  the  inhomogeneities.  This  procedure  can  be  used  to  randomize  the  travel  time  computation 
and  determine  fluctuations  in  signal  arrival  times  that  are  due  to  ocean  inhomogeneities. 

The  warming  of  the  upper  layer  of  the  sea  also  affects  acoustic  signal  travel  time.  For 
good  travel  time  resolution,  it  is  necessary  to  use  broadband,  high-frequency  signals.  To  avoid 
excessive  attenuation  due  to  absorption  and  surface  waves,  the  frequency  should  not  be  too  high. 
In  experiments  using  long  paths  in  the  ocean  (more  than  1000  km),  signals  of  several  tens  of 
hertz  are  used.  The  sound  signal  path  in  the  Fram  Strait  is  a  few  hundred  kilometers  long; 
therefore,  signals  with  frequencies  of  several  hundred  hertz  should  be  used.  The  sound  field  in 
this  frequency  range  normally  can  be  analyzed  by  ray  representation.  One  can  consider  signal 
travel  time  as  the  propagation  time  along  one  of  the  rays  between  the  transmitter  and  the  receiver. 
If  the  depth  of  the  transmitter  and  the  receiver  are  appropriately  chosen,  the  signal  travel  time 
along  one  or  several  rays  can  differ  significantly  from  the  travel  time  along  other  rays.  Thus,  if 
the  signal  frequency  band  is  relatively  broad,  it  is  possible  to  separate  the  signals  passing  along 
the  different  rays  by  delay  time. 
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A  numerical  ray-tracing  code  was  used  to  calculate  the  ray  structure  of  the  sound  field. 
This  provides  the  sound  field  level  through  both  coherent  and  noncoherent  summations  of  the 
sound  signals  propagating  along  different  rays  in  the  inhomogeneous  ocean.  The  code  takes  into 
account  spherical  and  cylindrical  spreading,  chemical  absorption,  surface  and  bottom  effects, 
the  sound  wave  attenuation  in  a  homogeneous  ocean,  and  anomalies  of  sound  signal  propagation 
due  to  inhomogeneities  of  the  medium.  At  program  input,  the  two-dimensional  sound  speed 
distribution  in  the  ocean  is  given,  and  the  positions  of  the  sound  source  and  receiver  are 
specified.  Propagation  times  for  the  sound  signal  along  the  pathway  were  calculated  for  different 
surface-layer  temperatures.  The  amplitudes  and  arrival  times  of  acoustic  signals,  propagating 
along  different  rays,  vary  as  the  environment  changes.  Stable  rays  sensitive  to  surface 
temperature  variation  were  chosen  for  the  simulations.  One  of  them  (ray  4s,  Fig.  7)  propagates 
initially  downward  and  experiences  four  reflections;  the  other  (ray  5s)  propagates  initially  upward 
and  has  five  contacts  with  the  sea  surface.  Three  receiver  positions  were  considered  at  depths  of 
100  (Ri),  200  (R2),  and  300  m  (R3),  respectively. 

The  ray  pattern  is  presented  in  Fig.  7.  The  rays  are  presented  for  signal  reception  at  levels 
of  100  and  300  m.  The  variation  of  the  hydrographic  environment  changes  the  ray  paths  and  can 
lead  to  the  disappearance  of  some  of  the  rays.  Ray  4s  at  receiver  R,  is  absent  for  the  initial 
environmental  condition;  it  appears  after  the  first  shift  and  disappears  again  on  the  fourth  shift. 
For  receiver  R3,  ray  4s  disappears  at  the  final  shift.  At  all  receivers,  ray  4s  disappears  at  the  four- 
step  shift  to  the  right,  while  at  the  three-step  shift,  this  ray  disappears  only  at  receivers  R,  and  Rj. 
In  all  the  other  cases,  this  ray  is  stable. 
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Figure  7.  Paths  of  two  rays  along  the  trace.  Ray  4s  contacted  the  sea  surface  four  times; 
ray  5s  contacted  the  sea  surface  five  times. 

Table  2  shows  the  mean  arrival  times  of  signals  propagated  along  ray  5s  for  receivers 
placed  at  depths  of  100,  200,  and  300  m,  along  with  root-mean-square  (rms)  values  of  arrival- 
time  fluctuations  averaged  over  1 1  realizations. 
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Table  2.  Arrival  Times  of  Signals 


Receiver  depth 

Ray 

5s  arrival  times 

rms 

(m) 

(s) 

(ms) 

100 

137.02741 

32.0 

200 

137.01450 

29.0 

300 

137.00314 

27.1 

Ray  5s  is  emitted  upward  from  the  source.  Table  2  demonstrates  that  arrival-time 
fluctuations  along  the  moderately  steep  rays  due  to  mesoscale  inhomogeneities  are  on  the  order  of 
30  ms.  These  fluctuations  are  easily  measured,  since  the  typical  rms  travel  time  uncertainty  for 
these  ranges  is  1-2  ms. 

Travel  times  for  sound  signals  propagating  along  different  rays  were  calculated  for  the 
different  temperatures  of  the  upper  200-m  subsurface  layer.  Typical  results  are  given  in  Fig.  8, 
where  the  amplitudes  and  arrival  times  of  the  signal  propagating  along  the  different  rays  are 
presented.  The  mean  arrival  time  variation  with  the  temperature  was  calculated.  It  was  obtained 
that  travel  time  variation  at  30  ms  (which  is  equal  to  fluctuation  rms)  corresponds  to  a 
temperature  change  at  0.6 °C. 

3.2  Sensitivity  of  Cumulative  Summation 

To  detect  the  temperature  trend  against  the  random  travel  time  fluctuations,  we  consider 
the  cumulative  sum  with  respect  to  the  total  train  of  arrival  times,  r,,  corresponding  to  the  signal 
propagation  along  the  different  rays.  When  one  arranges  the  signal  arrival  times  in  order  of 
increasing  value,  r,  <  /j  <  '3  ••  •  <  '/v^  the  cumulative  sum  can  be  written  as 


(3) 


The  normalization  factor  1  IN  is  used  to  compare  the  results  of  the  calculation  for  signal 
realizations  with  different  numbers  of  rays.  The  cumulative  sum  computation  procedure  is  as 
follows: 

•  The  signal  threshold  is  defined. 

•  The  appropriate  amplitude  of  the  signal  propagating  along  the  chosen  rays  is 
normalized  to  1  if  its  value  is  higher  than  the  threshold;  otherwise,  it  is  set  to  0. 

•  The  normalized  amplitude  of  each  successive  signal  is  added  to  the  sum  of  the 
normalized  amplitudes  of  the  preceding  signals  in  order  of  signal  arrival  time. 

•  The  total  sum  is  normalized  with  respect  to  the  number  of  terms  in  the  series. 
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Figure  8.  Signal  amplitude  and  arrival-time  dependance  on  the  subsurface-layer  temperature 
for  different  rays.  Distance  between  transmitter  and  receiver  is  200  km;  source  depth,  500  m; 
and  receiver  depth,  300  m. 


14 


The  following  analysis  of  the  cumulative  sum  is  based  on  the  computation  of  the  linear 
regression  coefficients  of  S(t)  with  respect  to  the  time,  t  ,  of  signal  arrival.  In  other  words,  the 
coefficients  a  and  b  are  defined  for  the  expression  at.  +  b,  which  gives  the  best  approximation  Sit.) 
in  accordance  with  the  least-squares  method. 


S  [s(t,) 


at 


(4) 


(  =  1 


This  procedure  is  outlined  schematically  in  Figs.  9a-c,  which  shows  the  arrival  signal  pattern 
Fig.  9a),  signals  that  exceed  the  threshold  after  amplification  and  limitation  (Fig.  9b),  and  the 
cumulative  sum  (solid  line.  Fig.  9c)  and  the  regression  line  (dotted  line.  Fig.  9c)  as  a  function  of 
time.  Regression  lines  corresponding  to  different  temperatures  of  the  layer  are  presented  in  Fig.  10 
as  a  function  of  signal  arrival  times.  In  Fig.  10,  the  labels  indicate  the  subsurface-layer  temperature 
increase,  AT  (AT  -  0  corresponds  to  the  initial  state  of  the  environment),  and  the  hydrographic 
pattern  shift.  So  (0  K,  0  km)  corresponds  to  the  initial  environment  realization,  (0  K,  30  km) 
corresponds  to  that  shifted  at  30  km,  and  (AT  K,  0  km)  corresponds  to  the  realization  with  increased 
subsurface-layer  temperature.  The  difference  in  slope  of  the  lines  corresponding  to  the  initial 
environment  realization  and  the  realization  that  shifted  at  30  km  indicates  the  influence  of 
background  noise  produced  by  mesoscale  inhomogeneities.  The  other  lines  correspond  to  different 
degrees  of  warming  of  the  upper  surface  layer.  It  can  be  seen  from  the  plot  that  the  slope  of  the  lines 
increases  as  the  temperature  of  the  ocean's  subsurface  layers  increases. 
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Figure  9.  Cumulative  sum  derivation,  (a)  Arrival  pattern;  (b)  signals  exceeding  the  threshold,  amplified 
and  limited;  (c)  cumulative  sum  (solid  line)  and  regression  line  (dotted  line)  as  a  function  of  time. 
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Figure  10.  Cumulative  sum  regression  lines  versus  signal  arrival  times  for  different  subsurface- 
layer  temperatures.  The  subsurface-layer  temperature  increase,  AT  (AT  =  0  corresponds  to  the 
initial  state  of  the  ocean),  and  the  hydrographic  pattern  shift  are  indicated. 


3.3  Strong  Mode-Coupling  Effect 

The  modal  approach  was  also  used  to  calculate  the  sound  signal  arrivals  (Shang  et  al., 
1996).  Numerical  simulations  have  been  conducted  for  the  following  four  cases: 


Case  1 
Case  2 
Case  3 
Case  4 


Lower  frequency  (/  =  20  Hz),  shallower  source  (Z,  =  50  m). 
Lower  frequency  (/  =  20  Hz),  deeper  source  (Z,  =  200  m). 
Higher  frequency  (/  =  200  Hz),  shallower  source  (Z,  =  50  m). 
Higher  frequency  (/  =  200  Hz),  deeper  source  (Z,  =  200  m). 


Results  are  presented  in  Table  3  and  in  Fig.  1 1 .  In  Table  3,  A^^  is  signal  amplitude;  t^  and  t 
are  the  arrival  times  of  the  signals  when  the  mode  coupling  effect  was  taken  into  account  and 
in  adiabatic  approximation,  respectively;  and  T  is  the  collective  arrival  time  of  the  signals 
(see  below). 


Ad 
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Table  3.  Modal  Travel  Times  for  Case  1  (/  =  20  Hz,  Z,  =  50  m)  at  224-km  Range 


t 

c 

T 

m 

m 

c 

Mode 

K 

(s) 

(s) 

(s) 

DC  =  0  m  s  ' 
(background) 


1 

0.80 

153.245 

153.243 

2 

1.0 

153.245 

153.217 

3 

0.37 

153.352 

153.302 

4 

0.42 

153.474 

153.327 

5 

0.43 

153.453 

153.401 

6 

0.98 

153.503 

153.472 

Collective:  "  153.309 

DC  =  4.5  ms-'  1  0.91  153.184  153.160 

2  0.71  153.145  153.122 

3  0.37  153.165  153.222 

4  0.62  153.349  153.252 

5  0.32  153.391  153.326 

6  1.0  153.453  153.419 
Collective:  153.287 

DC  =  9.0  ms-'  1  0.82  153.142  153.119 

2  0.66  153.088  153.072 

3  0.25  153.007  153.130 

4  0.59  153.275  153.185 

5  0.23  153.245  153.250 

6  1.0  153.399  153.358 
Collective:  153.224 

DC  =  13.5  ms-'  1  0.69  153.129  153.091 

2  0.57  153.074  153.046 

3  0.41  153.055  153.052 

4  0.40  153.244  153.119 

5  0.20  153.015  153.180 

6  1.0  153.336  153.290 
Collective:  153.180 
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Figure  11.  Modal  arrival-time  pattern  for  case  1  at  224-km  range. 
Mode  numbers  are  indicated  above  each  line. 
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3.4  Sensitivity  of  Collective  Arrival 

Strong  mode  coupling  takes  place  along  the  propagation  path  (see  Fig.  12).  For  this 
extremely  complex  propagation  environment,  the  arrival  times  for  a  fixed  mode  sometimes 
will  not  have  a  monotonic  relationship  with  the  surface  temperature  increase.  However,  certain 
kinds  of  averaged  arrival  times  could  have  a  monotonic  relationship  with  the  surface  temperature 
increase.  Here,  we  propose  the  collective,  or  amplitude-weighted,  arrival  time,  defined  by 


M 


^c- 

m    m 
m 

M 

m 
m 

(5) 


where  A^  is  the  modal  amplitude,  t^  is  the  arrival  time,  and  M  is  the  maximum  number  of 
modes  taken  into  account.  The  collective  arrival  time  for  case  1  is  listed  in  the  last  column  of 
Table  3.  Figure  13  illustrates  the  relationship  between  collective  arrival  time  and  temperature 
increase  for  all  four  cases. 

3.5  Current  Measurement 

To  determine  the  velocity  of  flow  across  the  acoustic  pathway,  the  scintillation  method 
can  be  used.  A  sound  wave  passing  through  the  ocean's  fine  structure  is  modulated,  producing 
an  irregular  pattern  of  amplitude  and  phase  perturbations  on  some  distant  receiving  plane.  These 
perturbations  are  produced  by  medium  changes.  The  evolution  of  the  scintillation  pattern  occurs 
because  of  advection  and  decay  of  the  ocean  inhomogeneities.  If  the  transit  time  of  the 
scintillation  pattern  across  the  detector  is  short  compared  with  eddy  lifetimes,  then  the 
intervening  transverse  flow  velocity  can  be  estimated  from  the  time  delay  of  the  correlation 
function  of  the  signals  propagating  across  the  flow  on  the  separate  paths  (Fig.  14;  Clifford  and 
Farmer,  1983;  Crawford  et  al.,  1990).  Further  investigations  are  needed  to  evaluate  the  feasibility 
of  this  method  in  application  to  a  large  distance  such  as  the  Fram  Strait  environment. 

Another  approach  of  transverse  flow  sounding  is  based  on  the  horizontal  refraction  effect 
(Voronovich  and  Shang,  1995).  Two  vertical  mode-resolving  arrays  should  be  used  to  receive 
signals  propagating  along  separate  paths.  The  phase  difference  of  the  modes  at  those  two  arrays 
is  directly  proportional  to  the  value  of  the  horizontal  refraction  caused  by  the  horizontal  gradient 
of  the  temperature  and  the  transverse  component  of  the  current  (Fig.  15).  Conventional 
tomographic  methods  can  determine  current  velocity  along  a  vertical  slice  between  two 
instruments.  When  combined  with  the  corresponding  temperature  measurement,  one  can 
determine  a  path-averaged  heat  flux  along  the  slice.  For  making  volume  and  heat  flux  estimates 
in  a  strait,  this  may  be  of  some  help,  but  the  method  requires  multiple  instruments,  essentially  to 
map  each  field  separately  in  a  tomography  mode.  One  possible  method  that  is  a  hybrid  point/ 
tomographic  method  would  be  to  place  multiple  inverted  echo  sounders  (MIES)  across  the  Strait 
in  a  zig-zag  pattern  to  measure  local  current,  temperature,  and  fluxes;  these  could  be  combined 
with  pressure  and  horizontal  electric  field  sensors  to  provide  additional  data. 
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Figure  13.  The  collective  arrival  time  versus  the  surface  temperature  increase  for  cases  1^. 
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Figure  14.  Ocean  current  measured  by  acoustic  scintillation. 
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Figure  15.  Ocean  current  measured  by  horizontal  refraction  tomography. 
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4.  PARAMETRIC  ARRAY  AS  A  MODE-SELECTIVE  SOURCE 

The  sound  propagation  in  a  complex  environment  such  as  the  Fram  Strait  is  characterized 
by  strong  mode  coupHng,  which  leads  to  mode  repopulation  (variation  of  the  modal  structure  of 
a  sound  field  with  distance).  The  arrival-time  pattern  varies  in  a  complicated  manner,  and  the 
arrival  times  for  a  fixed  mode  sometimes  will  not  have  a  monotonic  relationship  with  the 
subsurface-layer  temperature  increase.  The  application  of  a  parametric  array  (PA)  could  be 
useful  in  this  case.  A  parametric  array,  suggested  by  Westervelt  (1963),  is  based  on  the  nonlinear 
effects  of  sound  generation  as  a  result  of  interaction  at  a  high-frequency  wave.  Recent 
experiments  demonstrated  the  efficiency  of  a  PA  for  ocean  probing  (Andebura  et  al.,  1990).  The 
main  feature  of  the  PA  is  to  produce  a  narrow-beam,  broad-frequency-band,  sharply  directed 
radiation  (Andebura  et  al.,  1990)  that  could  be  used  for  selective  mode  excitation  in  a  stratified 
ocean.  Some  results  of  simulation  of  a  PA  application  for  the  Fram  Strait  sounding  are  presented 
in  the  following  figures.  In  Fig.  16,  the  modal  structure  of  the  sound  field  for  different 
orientation  of  the  PA  is  presented;  9  is  the  angle  between  the  axis  of  the  PA  and  the  horizon. 
The  sound  field  intensity  distribution  for  different  inclinations  of  the  PA  axis  in  comparison  to 
that  produced  by  a  point  source  are  given  in  Figs.  17-19.  It  is  seen  that  application  of  the  PA 
allows  sensing  of  different  depth  intervals  of  the  ocean.  In  Fig.  20,  the  collective  arrival  time  as 
a  function  of  the  ocean  temperature  increase  is  presented  for  a  point  source  and  for  the  PA  having 
different  orientations  (the  inclination  of  the  PA  axis  with  respect  to  the  horizon).  This  figure 
demonstrates  that  the  PA  provides  more  sensitive  probing  of  ocean  temperature  increases.  The 
PA  sound  field  features  could  also  be  used  for  current  measurements  (Dmitriev  and  Naugolnykh, 
1994). 


5.  CONCEPT  OF  STUDY 

5.1  Numerical  Modeling 

An  acoustic  monitoring  system  for  the  Fram  Strait  will  necessarily  be  based  on  an 
inversion  scheme  using  measured  travel  times,  intensity,  mode  structure,  and  refraction  angles  to 
obtain  temperature  and  the  current  field  and,  thereby,  the  heat  and  volume  flux  through  the  Strait. 
To  establish  the  inversion  scheme,  it  is  necessary  to  have  prior  information  available  about  the 
Strait  and  feasible  acoustic  propagation  models.  Inversion  schemes  are  built  upon  a  theoretical 
and  numerical  acoustic  propagation  model;  the  accuracy  of  the  inversion  result  depends  on  the 
accuracy  of  the  forward  model.  Therefore,  our  first  goal  is  to  do  computer  modeling  of  signal 
travel  time  variation  with  respect  to  water  temperature  and  current  using  the  acoustic  tomography 
and  scintillation  approach,  in  order  to  monitor  the  volume  and  heat  flux  through  the  Fram  Strait. 
Following  the  computer  modeling,  the  pilot  experiment  will  be  performed. 
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5.1.1  Physical  oceanographic  data 

Physical  oceanographic  data  would  best  be  made  available  as  a  small  database  on 
an  FTP  site.  Particularly  useful  data  for  our  planning  would  be  Atlantic/Arctic  waterfront 
characteristics,  mesoscale  eddy  field  characteristics,  internal  wave  field  characteristics,  and  ice 
cover  characteristics.  In  using  the  oceanography,  the  best  product  would  be  T,  S,  and  (m,v)  as 
functions  of  (x,  y,  z,  t),  where  /  includes  seasonal  effects.  Sound  speeds  are  easily  obtained  from 
these.  Also  useful  to  our  efforts  would  be  the  exact  bathymetry,  bottom  geoacoustic  properties, 
coastal  front  characteristics,  and  the  typical  wind  and  wave  fields. 


5.1.2  Forward  propagation  modeling 

First  consider  some  of  the  system's  characteristics.  Assume  that  the  system's  point 
sources  will  be  at  250  Hz,  with  a  100-Hz  bandwidth,  and  will  generate  192-dB  source  level 
initially.  (This  is  typical  of  HLF-5s,  which  were  used  very  successfully  in  the  Greenland  Sea.) 
Also,  assume  that  vertical  line  array  receivers  are  located  at  the  shallow  eastern  (Svalbard)  part  of 
the  system  (at  minimum)  and  probably  are  positioned  at  all  source  locations.  This  array  will 
nominally  span  the  water  column.  Further,  assume  that  computer  storage  of  data  at  the  receiver 
is  not  an  issue,  since  disk  storage  is  very  efficient  these  days.  Unless  we  use  cabled  systems, 
these  sources  and  receivers  will  have  limited  duration;  therefore,  temporal  sampling  strategies 
would  need  to  be  considered.  The  adequacy  of  the  signal  strength  received  and  the  stability  and 
identifiability  of  the  arrivals  will  also  need  to  be  considered.  Additional  points  to  address  are 
(1)  the  ice  scattering  and  surface  wave  rough  surface  effects,  (2)  bottom  interaction  use  of  bottom 
interacting  signals,  (3)  travel  time  perturbations  expected  due  to  temperature  effects, 
(4)  reciprocal  travel  time  signals  expected  due  to  current  effects,  (5)  scintillation  signals  expected 
due  to  advection  of  structure  through  the  tomographic  path,  and  (6)  horizontal  refraction  effects 
due  to  currents. 


5.1.3  Inverse  problem 

At  least  three  inverses  for  the  proposal  are  (1)  a  cumulative  inverse;  (2)  a  conventional 
tomography  inverse  for  T,  u,  v,  together  with  resolution  and  variance  estimates;  and  (3)  a 
rudimentary  scintillation  inverse  for  fluxes.  Additional  forward  modeling  and  inverse  schemes 
have  to  be  considered:  more  sophisticated  range-dependent  rays,  coupled  modes,  and  parabolic 
equation  calculations  for  forward  modeling,  as  well  as  horizontal  ray  refraction,  S-matrix,  and 
single-mode  generation  schemes  for  inversion  work. 


5.1.4  Observing  system  simulation  experiments 

In  the  long  term,  one  will  want  to  combine  the  results  of  measurement  with  ocean  general 
circulation  models  to  obtain  the  best  overall  estimate  of  the  four-dimensional  state  of  the  ocean 
(Evensen,  1994a,b).  The  data  would  consist  not  only  of  the  acoustic  observation  that  has  been 
the  primary  subject  of  this  report,  but  other  complementary  information  that  may  be  available, 
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such  as  the  results  of  satellite  observations  or  from  cable  voltages.  Observing  system  simulation 
experiments  (OSSEs)  must  be  performed  to  determine  what  is  the  best  mix  of  the  different 
measurement  techniques.  Some  of  the  possible  questions  that  need  to  be  addressed  are  (1)  what 
resolution  is  required  in  the  vicinity  of  the  different  fronts  and  at  other  locations;  (2)  how  well  is 
the  true  turbulent  heat  flux  measured  (if,  for  instance,  the  velocity  and  temperature  are  measured 
along  a  path  or  over  an  area);  (3)  how  many  degrees  of  freedom  are  necessary  to  describe  the 
flow;  (4)  what  are  the  correlation  length  and  time  scales;  etc. 


5.2  Pilot  Experiment 

5.2.1  Objectives 

The  specific  objectives  of  the  pilot  experiment  are  (1)  perform  scientific  studies 
on  the  feasibility  of  a  permanent  monitoring  of  the  Atlantic  water  and  heat  inflow  through  the 
Fram  Strait;  (2)  demonstrate  the  efficacy  of  a  combined  acoustics/conventional  sensor  system 
for  heat  flux  measurements;  (3)  obtain  a  full  year  of  oceanographic  data  to  initiate  the  monitoring 
effort;  (4)  test  some  novel  concepts  in  acoustic  monitoring,  which  may  allow  improvement  of  the 
system;  and  (5)  establish  requirements  for  needed  technical  equipment  for  long-term  monitoring. 

Proposed  tasks  are: 

Task  1.  Computer  modeling  of  signal  travel  time  variation  with  respect  to  water 
temperature  and  current  variations  using  archival  oceanographic  data  as  input. 

Task  2.  Develop  an  experimental  scheme  for  a  feasibility  demonstration  of  remote 
acoustic  monitoring  of  heat  and  mass  transfer  through  the  ocean  current  cross  section  based  on 
an  acoustic  scintillation  approach. 

Task  3.  Perform  a  pilot  experiment  to  measure  the  travel  time  variation  and  scintillation 
of  signals  crossing  the  Fram  Strait;  obtain  experimental  verification  of  the  principles  of  S-matrix 
heat  monitoring  (a  new  approach). 

Task  4.  Analyze  the  experimental  results  and  develop  recommendations  for  permanent 
remote  acoustic  monitoring  of  heat  and  mass  transfer  through  the  Fram  Strait  into  the  Arctic 
Ocean. 


5.2.2  Methods  of  investigation 

Different  methods  for  acoustic  monitoring  of  temperature  and  velocity  of  current  through 
the  Fram  Strait  will  be  used.  For  the  basic  acoustic  array,  three  modified  Avatoc  transceivers 
(Al,  A2,  and  A3)  and  one  DOAR  vertical-line-array  (VLA)  receiver  will  be  used.  The  array  is 
shown  in  Fig.  21.  The  two  eastern  receivers  will  allow  a  year-long  time  series  to  be  obtained  for 
(1)  standard  temperature  tomography;  (2)  reciprocal  current  measurements  (Al  to  A3,  A2  to  A3, 
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and  A 1  to  A2);  (3)  scintillation  measurements  (Al  to  A3/DOAR  and  A2  to  A3/D0AR);  and 
(4)  horizontal  refraction  measurements  (Al  to  A3/D0AR  and  A2  to  A3/DOAR).  Some  number 
of  conventional  oceanographic  sensors  are  needed  both  to  aid  the  acoustics  interpretation  and  to 
provide  data  for  the  overall  inverse.  To  help  look  at  temperature  near  the  fronts,  whose  east-west 
movement  across  our  acoustic  paths  is  critical,  T-strings,  inverted  echo  sounders,  and  satellite 
remote  sensing  will  be  used.  The  acoustic  Doppler  current  profilers  (ADCPs),  the  temperature 
and  conductivity  gage  (Seacats),  and  the  temperature  measuring  strings  (T-strings),  all  located  on 
the  acoustic  moorings  and  halfway  across  the  A 1 -to- A3  path,  help  provide  supplementary  flux 
measurements  to  verify  tomographic  measurements  and  help  with  the  acoustic  propagation 
analyses.  The  array  will  be  deployed  in  late  summer,  during  which  time  the  shipbome  sources 
and  other  instrumentation  to  test  the  various  flux-measuring,  inverse  schemes  proposed  will  be 
used.  This  deployment  cruise  would  last  2-3  weeks,  most  likely.  It  would  need  at  least  a 
medium-sized  vessel  accommodating  a  science  party  of  approximately  10-12  people.  The  basic 
array  of  sensors  would  be  left  in  the  water  over  the  next  year  and  recovered  the  next  summer. 
In  order  to  test  the  novel  inverse  schemes  in  a  timely  fashion,  we  would  have  to  transfer  data 
from  the  Al  and  A3  transceivers  and  the  DOAR  array  directly  after  the  pilot  experiment. 


5.2.3  Expected  results  and  the  project  period 

Expected  results  of  the  project  are  the  experimental  comparison  of  the  efficiency  of 
different  methods  of  acoustical  monitoring  in  variable  environments  typical  of  the  Fram  Strait; 
the  performance  of  a  feasibility  test  of  the  combined  acoustic-conventional  system  for  heat  flux 
measurements;  a  full  year  of  oceanographic  data  for  the  Fram  Strait  environment;  and  the 
original  data  to  design  an  optimal,  acoustical  measuring  system  for  sensing  of  volume  transport 
and  heat  flux  through  the  Fram  Strait. 


6.  CONCLUSION 

What  is  the  specific  problem?  Greenhouse  wanning  in  the  Arctic  is  expected  to  be 
two  to  four  times  greater  than  that  in  lower  latitudes.  If  this  happens,  it  could  significantly  raise 
global  sea  levels.  The  water  mass  and  heat  exchange  through  the  Fram  Strait  are  essential 
boundary  conditions  for  global  climate  models.  Previous  estimates  of  the  current  flux  of  the 
West  Spitsbergen  current  across  79°  N  vary  from  1.9  to  8  S^.  This  large  variation  in  water  mass 
exchange  introduces  errors  in  ocean  circulation  and  global  climate  models. 

What  are  we  doing?  An  acoustic  monitoring  scheme  is  being  designed  that  will 
transmit  low-frequency  sound  across  the  Strait  to  measure  average  temperature,  and  scintillation 
techniques  will  measure  the  transverse  current.  Combining  the  two  measurements  gives  an 
estimate  of  heat  flux.  Propagation  modeling,  already  under  way,  shows  the  acoustic  environment 
to  be  very  complex,  with  strong  mode  coupling.  A  parametric  array  could  be  used  to  obtain  fine 
depth  resolution,  and  the  low  acoustic  power  in  a  required  low-frequency  range  will  make  the 
measurement  "whale  safe." 
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What  is  the  payoff?  Improved  measurements  of  Arctic  warming  will  reduce  the 
uncertainty  in  the  global  warming  models.  Data  assimilation  taken  during  the  Fram  Strait 
monitoring  will  provide  an  essential  contribution  to  the  ATOC-Arctic  project  and  will  contribute 
to  a  reduction  in  the  uncertainty  in  the  modeling  and  prediction  of  global  climate  warming. 
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